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EMERIES EROM EARTH ORBIT THAT SIMULATE ACCELERATION 


STRESSES OF MANNER PLANETARY MISSIONS 
By Frederick C. Edwards 
Ames Research Center 


SUMMARY 


The ability to duplicate the acceleration stress during a portion of the 
earth return of a planetary mission vehicle by appropriate control of an earth 
orbital vehicle is investigated. The results show that during the entry of 
an orbital spacecraft it is possible to approximately match the acceleration 
profile corresponding to the initial dive, the capture maneuver, and part of 
a constant altitude deceleration phase of a planetary mission return entry. 

A two-step approach to the solution is examined: The analytical step results 
in a closed-form approximate solution for the lift control. An analog com- 
puter simulation of the entries uses results from the analytical step and 
gives a more precise solution for the control. It also permits an evaluation 
of the cost to the earth entry vehicle in terms of the aerodynamic heating, 
range capability, and retro velocity requirements. 

Emphasis is placed on simulating the planetary mission entry at twice 
orbital velocity, but the results also apply to entry at other supercircular 
velocities. An Apollo-type lifting entry vehicle was used for the example 
entries presented. 


INTRODUCTION 


The long duration missions to the planets would require astronauts to 
endure far longer periods of weightlessness (300 to 600 days) and higher 
reentry acceleration levels (lO g and greater) than has been experienced in 
past earth orbital missions. If the present man-system relationship is main- 
tained, the astronaut will be required to monitor, guide, and control during 
entry from these missions. The effects of these extended missions and, in 
particular, the weightless environment on his ability to tolerate high entry 
stresses and at the same time perform control functions must be investigated 
prior to the planetary missions. 

Entry from long duration earth orbit missions provides opportunities for 
studying performance of astronauts after they have been exposed to extended 
periods in a weightless environment. This report indicates a means whereby 
earth orbital spacecraft may be used during the entry phase to give meaningful 
data on the effects of the mission environment on the performance of the 
astronaut during planetary entry; the technique is to control an earth orbital 
vehicle during entry so as to duplicate a planetary mission entry 
acceleration -time profile. 


The investigation initially considers the flight dynamics of a planetary 
mission return vehicle (to be called the reference vehicle) , followed by a 
discussion of the general character of the entry of an earth orbital vehicle 
(to be called the flight vehicle) when it is controlled to duplicate the plan- 
etary mission entry acceleration profile. An approximate equation is 
developed for the required lift control. Numerical results are obtained from 
the approximate equation and from a more complete analog computer simulation 
of the trajectory motions. The effects of acceleration control on the vehicle 
range capability, aerodynamic heating, and retro velocity increment are deter- 
mined. A roll -modulated lifting capsule (Apollo) is the spacecraft considered 
in this study. 

SYMBOLS 

At total acceleration, g 

D 

C D drag coefficient, y p yg s 

C L coefficient, y p - g£ , 

D drag force , N 

g local value of gravity acceleration, m/sec 2 

h altitude, m 

H total integrated heat load, j/m 2 

Ki ,K 2 control equation gain constants 

L lift force, N 

m mass of vehicle, kg 

q heating rate per unit area, W/m 2 

r distance from planet center, m 

S reference area for drag and lift coefficients, m 2 

S D dov/nrange distance traveled, km 

crossrange distance traveled, km 
t time, sec 

u circumferential velocity component normal to radius vector, m/sec 

Ue circular orbital velocity \[ gr , m/sec 


2 



u 

V 

V 
AV 

^1 

w 


dimensionless velocity ratio ^7 
resultant velocity v u 2 + w 2 , m/sec 

x. 

Uc 

dimensionless change from initial velocity. 


Vi 


V 


s/gr 


Xh^h.Zh 
Z 


dimensionless velocity at peak acceleration 
vertical velocity component (along direction of radius vector), m/sec 
physical components of external force, Ef 


P 

7 

■n 

p 

T 

9 


_ Ct)S ct 

dimensionless function of V given by pV 


atmospheric scale height, l/m 

flight -path angle given by arc tan 

retro thrust vector direction, deg 

atmospheric density, kg/m 3 

dimensionless time, /f— t 

7 v x o 

roll angle, deg 

lateral deflection angle, deg 


—w 

u 


Subscripts 

com command 
nom nominal 

ref reference (time function) 

0 reference (constant value) 

1 flight vehicle (initial velocity = l) 

2 reference vehicle (initial velocity = 2) 
f final value 

i initial value 
I value at maximum acceleration 
v vert ical component s 
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I 


hi 1 1 ii 


i ii ■ 



C) 

0 ) 

n 


Superscripts 


differentiation with respect to time 
differentiation with respect to velocity 
dimensionless quantity 


ANALYSIS 


General Remarks 


This section will discuss the control required, during the atmosphere 
entry of an earth orbital mission vehicle, to duplicate the acceleration time 
history during the atmosphere entry of a planetary mission vehicle. A two- 
step approach is presented. The initial step considers only one component of 
the entry vehicle motion and leads to an approximate but useful analytical 
solution. A second step, employing an analog simulation which includes three 
degrees of vehicle translational motion and one degree of body rotational 
motion, uses the approximate results from the first step to obtain a more 
precise solution. 


Planetary return 


mission return (fig. 1 


Orbital 

return 


Planetary 

return 



(a) ph vs. V 

Figure 1.- Typical earth entry trajectories for 
planetary and earth orbital mission vehicles. 


reference ve hicle) . - The trajectory for a planetary 
a^n is assumed to have an entry velocity twice the 

local satellite velocity. The par- 
ticular entry shown is one of sev- 
eral types that could be considered 
and is referred to as a constant 
altitude entry. To achieve a suc- 
cessful entry the vehicle must enter 
within the safe corridor between the 
undershoot and overshoot boundaries 
of figure l(a) . During the initial 
dive into the atmosphere the vehicle 
acceleration increases. At maximum 
acceleration (pull out) the entry 
control system (either manned or 
automatic) is required to perform a 
"capture” maneuver to insure that 
the vehicle does not skip out of the 
atmosphere or dive in to such a 
degree that destructive acceleration 
forces build up. In figure l(a) 
the vehicle enters near the under- 
shoot boundary. It is stressed 
in reference 1 that entry near this 
boundary presents a critical con- 
trol problem. The vehicle lift 
vector will be initially full up to 
insure that the maximum acceleration 
limit is not exceeded, then, 


k 



at peak g, the vehicle is rolled over so that lift and gravity forces balance 
the centrifugal force. The next portion of the entry is a constant altitude 
deceleration phase until sufficient lift can no longer be generated, at which 
time an equilibrium glide with lift vector full up is flown to landing. 

Orbital return (flight vehicle) .- The entry trajectory for return from an 
earth orbital mission (fig. 1(a)") is the particular one that duplicates the 
acceleration time history for the planetary mission return vehicle. A 
comparison of the altitude velocity profiles indicates how this match is 
achieved. 


The orbital return vehicle dives to a lower altitude (higher atmospheric 
density) than the planetary mission return vehicle, and after pull out, pulls 
down into the atmosphere at an increasingly greater rate as the velocity 
decays. At each instant when the accelerations along the trajectories are 
equal, the velocity of the planetary return vehicle will exceed that of the 
orbital return vehicle. It should be clear that in order for the aerodynamic 
accelerations to be equal for the two vehicles, the effect of the higher 


velocity of the planetary return vehicle 



(b) J3h vs. AV 
Figure 1.- Concluded. 


is being compensated by flying the 
orbital vehicle at a lower altitude 
where the atmospheric density is 
greater. With this technique, the 
dynamic pressures [(l/2)pV 2 ] for the 
two vehicles are equivalent. A more 
direct comparison of the two entry 
trajectories is shown in figure l(b), 
where the two entry trajectories 
presented in figure l(a) are now 
plotted_as a function of a new vari- 
able AV, the total velocity change. 1 
(in this plot the trajectories pro- 
ceed toward the right whereas in 
figure l(a) they proceeded toward 
the left.) In figure l(a), the two 
trajectories started at different 
velocities ( (Vi ) i - 1-0, (V 2 )j_ =■ 2.0) 
and proceeded to zero velocity at 
the terminal point. In figure 3_(b), 
they start at the same point, AV - 0, 
and_end at different points 
((AV x ) f - 1.0, (AV 2 ) f - 2.0). The 
difference in altitude between the 
two entry trajectories is more 
readily apparent in this figure than 
^-The difference between the 
initial velocity of the vehicle and 
the velocity at appoint along the 
trajectory (AV - V± - V) will be 
equivalent at equal times for both 
vehicles if the acceleration profiles 
are matched. A further discussion 
of AV will appear later in the text. 
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in figure l(a). In the lower part of figure l(b), traces of the total accel- 
eration for the two entry trajectories_are presented. Note that the accelera- 
tions are equal during the range of AV from 0 < AV < 1 . 0 , causing the two 
plots to appear as a single trace. It is during this increment in AV that 
it is feasible to match the acceleration profiles. 


Assumptions 

The entry vehicle assumed for this study was a blunt -faced lifting cap- 
sule trimmed at the appropriate angle of attack to give a desired L/D. This 
trim condition is set prior to entry by offsetting the center of mass from 
the center line of the vehicle. During the actual flight there is no direct 
control of the magnitude of the lift. The direction of the lift vector is 
controlled by rolling the vehicle about the velocity vector, giving rise to 
components of lift, and thus control, in the vertical and lateral directions. 


It was assumed that the vehicle aerodynamic coefficients did not vary 
with Mach number and Reynolds number and that the L/D and ballistic coeffi- 
cients were only functions of trim angle of attack. These aerodynamic coeffi- 
cients are given in table I as a function of trim angle of attack. 


The vehicle is assumed to enter the atmosphere of a nonrotating cylindri- 
cal earth under the influence of aerodynamic and gravitational forces. (An 
exponential approximation for the 1962 ICAO standard atmosphere is used 
throughout.) The differential equations describing the vehicle motions during 

entry are found by summing the accel- 
erations in the vertical, horizontal, 
and lateral directions at an arbitrary 
point along the trajectory. The fol- 
lowing nonlinear differential equa- 
tions have been derived in refer - 
2 and 3 • 


TABLE I.- AERODYNAMIC PARAMETERS FOR SIMULATED 
VEHICLE. 


a, deq 

L/D 

m/C D S, kg/m‘ 

0 

0 

269 

6 

0.10 

273 

15 

.25 

293 

31 

.50 

391 

47 

.70 

742 


ences 


d 2 h 


V 2 


cos 7 



dt 


dt 2 


D 

m 


sxn 7 - Ur 


cos 7 


(1) 


V cos 7 


D 

tn 


dh 

dt 


=- ) sin 7 + cos 7 


d 2 S^ 


dt 2 m 


L 

=• — cos 
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( 2 ) 


( 3 ) 
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In this study, as noted previously, solutions are obtained in two steps. 
Initially, only motion in the vertical direction is considered and equa- 
tions ( 2 ) and ( 3 ) are disregarded. An approximate but useful solution is 
obtained for the vehicle (L/D)v time profile. The second step of the 
approach considers all three equations of motion. These are programmed in a 
modified form and solved using an analog computer simulation. In addition to 
the (L/D) v time profile, more extensive results are obtained from this step 
of the approach. 


Control Technique 

In this study, an orbital mission return vehicle is controlled in a 
closed-loop manner during entry to match its flight trajectory variables to 
the corresponding reference trajectory variables. The reference trajectory 
variables are those generated during the entry of a planetary mission return 
vehicle. The variables controlled are the acceleration and the rate of change 
of acceleration. The control equation 

cp(t) = q>(t) + KiUtOt) - A t (t)] + K 2 [i t (t) - j t (t)] (10 

com nom fit ref fit ref 

is solved continuously during the flight trajectory to determine the value of 

the control cp(t) required to match the acceleration along the flight and 
com 

reference trajectories. The terms Ki and K 2 are constant feedback gains in 
equation (4). The term cp(t) is the nominal roll -angle profile along the 

nom 

flight trajectory. This is not equal to the nominal roll -angle profile of the 
reference trajectory because the initial velocity of the flight trajectory is 
only half the velocity of the reference trajectory. Although the acceleration 
profiles for the two trajectories can be Identical over a portion of the entry, 
the control profiles required to achieve this identity will be quite different 
because of the difference in velocity that exists throughout the portion of 
the entry where an acceleration match Is achieved. Consequently, to proceed 
with the analog solution, the nominal control profile for the flight trajec- 
tory must be determined with sufficient accuracy so that the linearized 
control equation is valid. 

Since the reference trajectory nominal profile ((L/D)^y versus velocity) 
and the initial velocity for the two trajectories are known, the corresponding 
flight trajectory control profile can be determined. It will be shown that 
an analytical approach results in a closed-form approximate solution for the 
flight vehicle control profile and allows the use of the linearized control 
equation in the analog computer approach to finding a more precise solution 
to the control profile. 


Analytic Approach 

To determine the proper trajectory control, during the entry of an earth 
orbit mission vehicle, to match the acceleration profile of a planetary entry 
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vehicle we shall initially use an approximate hut direct method of evaluating 
the control. The two entry vehicles may have dissimilar aerodynamic and mass 
characteristics but the total accelerations at an arbitrary instant of time 
must be equal. Thus, to indicate the variables that influence the match: 

At! 



where subscript 1 refers to the vehicle entering at local satellite velocity 
and subscript 2 refers to a planetary return entry at twice local satellite 
velocity- Equation (5) ma y be expressed as the difference in altitude between 
the two entry trajectories at each instant of time (which is necessary in 
order that their accelerations be equal) by assuming an exponential variation 
of density with altitude (p - Po^"^* 1 ) . The resulting expression is 


k-2 





(CpS/mg) g / l + ( L /D)/ 
(CpS/mg ) 1 J 1 + ( L / D ) i 2 


(6) 


Equation (6) indicates how the altitude profile of the two entry trajectories 
should differ as a function of the velocity ratio V 2 /Vi and the aerodynamic 
and mass characteristics of the vehicles. The second term on the right is 
constant during entry; thus, the altitude difference varies exponentially 
with velocity ratio only. This exponential variation is apparent in fig- 
ure l(b) for that portion of the orbital return trajectory (after pull out) 
corresponding to the constant altitude portion of the planetary entry 
trajectory. 

It may be noted that equation (6) could be a solution to the differential 
equation of motion (eq. (l)), except that the independent variable is velocity 
instead of time (as in eq. (l)), and the dependent variable is the difference 
in altitude rather than altitude directly. If the two equations were com- 
patible, the conditions under which equation (6) would be a solution to the 
differential equation (l) could be determined by directly substituting the 
solution into the differential equation and determining the values of (L/D)^ 
required to solve the resulting equation. Assuming for the moment that in 
equation (6) the reference altitude profile (h 2 vs. V 2 ) is known, we shall 
check to see for what conditions a solution is obtained. First, the indepen- 
dent variable must be changed so that the differential equation and its solu- 
tion correspond. The independent variable chosen will be neither the time nor 
velocity, as in equations (l) and (6), but the new variable (designated the 
dimensionless total velocity change) defined as 
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(7) 


AV = Vi - V 

where Vj_ is the initial velocity_of the vehicle, V, the velocity at any time 
along the entry trajectory, and AV, the difference between these two values. 
It is proposed that the values of this new variable will be common to both the 
reference and flight trajectories, in other words, the total changes in veloc- 
ity at a corresponding time along the trajectories will be equal if the total 
acceleration profiles are precisely matched. This requires that the integrals 
of the two total acceleration profiles be equal: 

^ = C Afci V = Jj \ s *r’ (8) 


Actually, equation (8) is only approximately true. The conditions under which 
it is used are presented in appendix A. 

Introducing the new variable AV into equation (l) and making the small- 
angle approximation cos y=l and sin 7=7 results in the simplified 
expression 


_ d(Vy ) = 1 - V s 

dAV A t /s/l + (L/D) 2 



(9) 


where the rate of climb, V7, has been substituted for its equivalent, — . 

v/gr at 

The control variable to be evaluated, (l/D) v , is the vertical component of 
the total lift-to-drag ratio of the vehicle. 

In equation (9) the derivative term is evaluated by differentiating 
equation (6),_a solution to the differential equation. Differentiating with 
respect to AV and substituting in terms of the rate of climb variable Vi7x 
yields an expression for the flight -path angle 73.: 


Vi7i 


V 2 7 2 


ll + (L/D)p 2 
1 + (L/D)x 2 


2At 

prVxV 2 


( 10 ) 


After a second differentiation. 


<a.Vx7i 

_ av 272 

/ 1 + (l/d)/' 


’l_ 

1 

PrV 2 7 2 

dAV 

dAV N 

/ 1 + (L/D)x 2 

prVxV 2 

Ji 

" v 2 

At/s/l + (L/D) 2 2 J 


where subscripts 1 and 2 correspond to the values. for the flight and reference 
vehicle, respectively. The total acceleration A^. has no subscript since it 
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is equivalent for both vehicles. Equations (lO) and (ll) axe substituted 
directly into the differential equation (9)> which is solved for the required 
value of (l/D)iv which makes equation (6) a solution of equation (9)- The 
result is: 


/l\ 1 - v x 2 | A + (l/d) 2 2 

A [A 1 - v 2 2 

At /A + (L/D)i 2 J 1 + (L/D)! 2 

k t /sfl+ (L/D) 2 2 


+ ZB- 

Vi 


'i + (l/d) 2 2 

1 + (L/D )i 2 



Vi + 2 


n/i + (L/D)! 2 ]- 


2K t (2V 2 - Vj) 

Pr(ViV 2 ) 2 ~ 


( 12 ) 


Substituting for the velocities _Vi and V 2 in_terms of the common variable 
AV and the Initial velocities (Vx)i = 1 and (V 2 )i = 2 gives the expression 
for the vehicle control: 


L\ _ 2AV - AV 2 
' D Av At/s/l + (L/D)! 2 



Uav - av 2 - 3 72(1 + 2J1 + (l/d ) 2 2 

A t //l + (L/D) 2 2 (1 - AV) 


2^(3 - AV) 
pr[(l - AV) (2 - AV)] 2 


(13) 


In equation (13) the lift required for control of the flight vehicle is 
related to the lift characteristics of both the flight and reference vehicle 
and the trajectory parameters of the reference vehicle only. The control is 
a direct function of the vertical component of the L/D of the reference 
vehicle , (l/D) 2 v, the L/D of the reference and flight vehicles, the integral 
of the total acceleration AV, the magnitude of the acceleration, At, and an 
indirect function of the rate of change of acceleration through the flight - 
path angle, Although it is not apparent at this point in the discussion, 

all of these parameters are known quantities. Thus, we have determined that 
equation ( 6 ) will be a solution to the differential equation if the vertical 
component of the lift vector varies as specified by equation (13) • Equa- 
tion (13) is then the closed-form approximate solution for the flight vehicle 
control problem. In order to solve equation (13) quantitatively, thejrefer- 
ence trajectory must be specified in terms of the variables (l/d) 2 vj At, and 
72 as a function of the independent variable AV. It shall now be indicated 
how these quantities may be determined and a numerical solution obtained. 
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The altitude versus AV profile for the reference trajectory (fig. 1(b)) 
is conveniently divided into three distinct segments. The initial dive into 
the atmosphere up to the point of pull out is the first phase. The constant 
altitude portion is the second phase, followed by the third, designated an 
equilibrium glide. Our interest will be in duplicating the acceleration pro- 
file of the first segment and as much as possible of the second phase. The 
reference variables for these two segments will be specified separately. 

Initial dive phase . - The initial dive into the atmosphere may be consid- 
ered that portion of the trajectory from atmospheric encounter (i.e., 

Af- = 0.05 g) until maximum acceleration. This portion of the trajectory has 
been solved in closed form for restricted cases by Gazley, Loh, Lees, and 
others and, in a more general form, by Chapman. Chapman (ref. 2) reduces the 
two basic differential equations for longitudinal motions (eqs. (l) and (2)) 
to a single approximate equation. This single equation is a nonlinear second- 
order differential equation in terms of a generalized dependent variable Z 
and the independent variable velocity, u. Although Chapman’s equation 
requires numerical machine calculations, the general results obtained apply to 
problems involving a range of body shapes, mg/CpS, and planetary atmosphere,". 

A Set of precomputed entry trajectories, presented in terms of generalized Z 
functions as well as the more useful values of At, 7, and V, is presented in 
tabular form in reference 4 for several different initial velocities and L/D 
values. The values from these tables may be used successfully as reference 
trajectory variables for the initial dive into the atmosphere. The accelera- 
tion, flight -path angle, and velocity for an initial velocity of V± - 2.0 
and a specified (l/d) 2V , (l/D) 2 , and (L/D)! may be taken from the appropriate 
table_of reference 4 andjsubstituted directly into equation (l3) as a function 
of /$f (where AV = 2 - V) . Equation ( 13 ) may then be solved for the (l/d) 1v 
of the flight vehicle and plotted as a function of velocity or the correspond- 
ing time (also given in the tables) . 

Although the above procedure is feasible, it was not directly followed 
for the example presented in figure 3; a new table of values was computed 
(from Chapman’s equations) that resulted in a trajectory with a peak accelera- 
tion of precisely 10.0 g. This figure will be discussed in a later section. 


Constant altitude phase . - For constant altitude flight the computation of 
the required (L/D ) 1 is jgreatly simplified because both the reference flight - 
path angle 72 and d 2 (V 7 )/dAV^ are zero. Thus the required (l/D) 2 v> evalu- 
ated from equation ( 9)5 is 



1 - v 2 2 


VVl + ( L /D); 


4av - aV* 2 - 3 

A t /s/l + (L/D) 2 2 


(14) 


When equation (l4) is substituted into equation ( 13 ), the terms in the bracket 
will be identically zero and the equation will reduce to the simpler form: 



2AV - AV 2 
AtA/l + (L/D K 


2At(3 - AV) 

0r[ 1 - AV) (2 - AV)] 2 


(15) 
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In equation (15) only the reference acceleration profile need be speci- 
fied, which, for the constant altitude (p = const) segment, will vary as the 
square of the velocity: __ 

A - s h(0 * (2 - 

where Atp/Vi 2 is a constant and is evaluated at pull out (maximum accelera- 
tion) . Substituting equation (l 6) into (15) and simplifying yields 



Vj Jl + (l/D)! 2 / AV \ A-fcj 

I 

co 

Atj \2 - AV/ 2pr 

L(i - a\0 2 J 


(17) 


The solution of equation (17) for the constant altitude portion along with 
(13) for the initial dive, comprises the approximate solution to the problem 
of matching the acceleration profiles of the two vehicles, one entering at 
planetary mission return velocity (Vj_ = 2), the other at orbital return veloc- 
ity (Vj_ - 1). From this approximate solution a precise solution can be 
obtained on the analog computer. 

The problem discussed earlier in specifying the nominal roll profile 

cp(t) in the vehicle control equation (eq. (4)) may now be resolved since a 
nom 

determination of the (L/D) 1v profile is, in effect, a determination of the 
roll profile (since cp = cos " x [ (l/D) v /(l/d) ] • An inconsistency exists here 
though, since the solution of equations (13) and (17) gives (l/D)iv as a 
function of the velocity variable AY and not of the variable time. This is 
a problem for only the constant altitude portion of the trajectory since the 
tabulation of the reference trajectory variables in reference 4 for the ini- 
tial dive includes the variable time. For the second phase of the entry 
trajectory a relatively simple relationship between time and AV is derived 
in appendix B and is presented here. 

t ‘ r ji (s-h - rhfi) * *1 (l8) 

The quantities Atj, AVj, and tj are the acceleration, velocity change, and 
time, respectively, at initiation of the second phase (i.e., pull out). 

These values are numerically equivalent to those existing at the end of the 
initial dive. Using equation (l8) in conjunction with equation (17) permits 
the evaluation of the nominal roll profile as a function of time. The profile 
in this form may be used in the control equation (eq. (4)) as the nominal 
roll profile. 


Analog Computer Approach 

The analytical step of the approach presented in the previous section 
gave an approximate solution for the L/D versus time (control) profile. To 
obtain this solution, it was assumed that the flight -path angle was small 
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and that the three vehicle translational equations were independent so that a 
solution based on a single differential equation of motion could be obtained. 
We shall now proceed with the more general analog computer step of the 
approach in which the above assumptions are not necessary. 

For this step, the full set of the three differential equations was sim- 
ulated on the computer. This allows a more precise solution for the control 
profile, provides flexibility for generating the desired reference trajec- 
tories, and permits evaluating the costs of performing the acceleration match 
(in terms of the sacrifice in vehicle ranging capability and changes in the 
aerodynamic heating on the vehicle) . The set of differential equations in the 
form presented in the previous section (eqs. (l), (2), and (3)) are not suit- 
able for direct programming on an analog computer. The large range of the 
program variables , plus the fact that the net radial acceleration is at times 
a small difference of large forces, presents scaling problems on the analog 
computer which greatly affect the accuracy of the results. However, if a more 
convenient axis system is introduced and the equation is modified to a pertur- 
bation form, the difficulties are reduced. The use of a modified flight -path 
axis system, introduced by Fogarty and Howe (ref. 3)j will allow continuous 
simulation from orbital velocity and altitude through entry without rescaling. 
The modified equations in this form are presented in appendix C. The set of 
equations represents the three degrees of vehicle translational motion to 
which is added one degree of body rotational motion (roll axis) and additional 
auxiliary relationships to give the aerodynamic heating and to solve the 
vehicle control equation. 


Circular velocity entry 
Initial conditions 


Vi n h t 



Vi r L 


Supercircular velocity entry 
Initial conditions 


Figure 2.- Block diagram of analog computer 
study . 


The analog computer set-up is 
represented diagrammatically in fig- 
ure 2. The two large blocks repre- 
sent the respective trajectory 
dynamics for the two vehicles. Two 
identical sets of the translational 
differential equations of motions 
were programmed on the analog com- 
puter. One set, representing the 
motions of the planetary mission 
entry vehicle, generates the refer- 
ence trajectory state variables and 
removes the need for programming a 
precomputed reference. The second 
set represents the motion of the 
orbital mission entry vehicle. The 
two sets of equations are solved 
concurrently; thus, at each instant 
of time, the corresponding values of 
the state variables are available 
for comparison to give the error 
signal. The trajectory of the ref- 
erence vehicle is completely indepen- 
dent, that is, different reference 
trajectory profiles are found through 
proper choice of the initial 





conditions. The trajectory motions of the flight vehicle are related to those 
of the reference vehicle through the control equation (eq. (4)). The mechani- 
zation of the control equation is also shown in figure 2. The state variables 
generated by the trajectory dynamics are compared to obtain the error quan- 
tities (Mt and M^.) ; these are weighted appropriately (empirically) and 
summed with the nominal roll profile to obtain the commanded roll angle. The 
flight vehicle is commanded to roll so as to reduce the error quantities. 

The approximate solution determined from the analytic step of the 
approach of the previous section was used as the initial estimate of the nom- 
inal control in the mechanized control equation. Since this cgntrol profile 
is only approximate, an error in the control parameters (Mt* would 

develop even if no trajectory disturbances occurred. For the computer mecha- 
nization, provisions were available for manually adjusting the shape of the 
approximate nominal profile. By proper adjustment so that the magnitudes of 
the error quantities were reduced, a better approximation of the solution for 
the nominal control profile could be obtained. If the error quantities could 
be reduced to zero, then the true nominal and, consequently, the exact solu- 
tion for the control would be found. The solution was not pursued to this 
extent since the lift required for certain regions of the control profile 
could not be achieved during the entry of an actual vehicle. Therefore, a 
compromise of the original objective of obtaining a precise match was required. 
The nature of this compromise will be explained fully in the next section. 


DISCUSSION AND RESULTS 


General Remarks 

In the analysis section, an analytical expression (eq. (13)) was derived 
which gives an approximate solution for the control of a vehicle entering 
from orbit so that the acceleration profile that results would match that of 
a planetary mission vehicle entering the earth 1 s atmosphere at twice orbital 
velocity. In this section, it will be shown that although It is Impractical 
to obtain a precise match of the acceleration profile, an approximate match 
is obtainable and has certain advantages which will be expounded. This will 
be followed by an extensive discussion of the relationship that must exist 
between the orbital vehicle aerodynamics and trajectory characteristics and 
those of the planetary entry vehicle In order to obtain an approximate match. 
Finally, the cost of performing the acceleration match in terms of the result- 
ing aerodynamic heating on the vehicle, the constraints on the vehicle maneu- 
verability for range control, and the changes in the required velocity 
increment for the retro maneuver will be discussed. 


Matching Considerations 

In figure 3; a typical set of entry trajectories is presented for which 
the flight vehicle is controlled to generate the acceleration profile of the 
reference vehicle. The altitude, acceleration, and_lift control profiles are 
presented as a function of the velocity variable AV. The reference 




Planetary mission 

entry vehicle (reference) 

Earth orbital mission 

entry vehicle (flight) 



Figure 3.- Typical variation of earth entry 
trajectory parameters with velocity during 
planetary mission entry simulation. 


acceleration profile has a peak of 
10 g and could correspond to the 
entry of a planetary mission vehicle 
near the 10 g undershoot boundary. 
Figure 3 shows that the reference 
vehicle lift vector is up (cp = 0°) 
during the initial dive (dashed 
curve) , rolled over near peak acceler 
ation to a downward direction to 
insure capture, and modulated for 
equilibrium during the subsequent 
period of constant altitude 
deceleration. 

Two entries of the flight 
vehicle are illustrated. One entry 
is presented in which a precise accel 
eration match is achieved. This 
means that the acceleration is equal 
to the reference value at each point 
along the trajectory. Thus, the 
acceleration trace is identical to 
the reference acceleration trace 
(dash curve). The second entry repre 
sents an approximate match, which 
means that during the initial dive . 
only the maximum values of A^ and At 
are matched rather than the accelera- 
tion at each point. During the sub- 
sequent deceleration, the flight 
vehicle is controlled to achieve a 
close match of the reference accelera 
tion profile in a closed-loop manner. 


In figure 3 , the profiles for the reference trajectory were generated by 
programming Chapmans entry equations on a digital computer. The control pro- 
file for the precise match is_ the quantitative solution of equation ( 13 ) using 
the outputs of At, 72 , and V 2 from the digital program as the reference. 

The altitude profile is computed using equation (6) of the analysis section. 
Figure 3 shows that the L/D required for a precise match as computed from 
equation ( 13 ) is too large, at both the onset and end of the entry, to be 
achieved by the vehicle considered here. Although a precise match of the plan 
etary mission entry acceleration profile by an earth orbital mission vehicle 
is physically unrealizable with a capsule -type configuration, because of the 
excessive aerodynamic lift required, the acceleration profile could be matched 
over portions of the entry. Figure 3 shows an example of this approximate 
match. The profiles for the approximate match were generated on an analog 
computer simulation of the entry with the same reference acceleration profile 
as that used to obtain the control profile for the precise match just 
described. Thus, this reference is also the dashed curve (A^ vs. AV) on fig- 
ure 3* The approximate solution has the characteristics of matching the mag- 
nitudes of the peak acceleration and the peak acceleration rate during the 
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initial dive phase of the entry while the control profile retains the general 
shape of the control profile corresponding to the planetary mission entry tra- 
jectory. The approximate match is achieved with constant L/d (constant roll 
angle) during the initial dive, a roll-over maneuver at peak g, and then 
closed-loop control of the lift vector about the nominal during the subse- 
quent period of acceleration. During the terminal portion of the entries 
(AV >0.5) the magnitude of the roll angle progresses toward 9°° (zero L/D) 
on fig. 3(b)) for the planetary mission entry, while for the orbital mission 
vehicle, the roll angle progresses toward l80° (negative L/d) as shown in 
figure 3(b). The difference in the two acceleration profiles in figure 3(a) 
is deceiving. Actually, the match on a time history plot is much closer than 
it appears here on a velocity history. The increasing error apparent in the 
two acceleration profile^ as AV iuiereases reflects the inaccuracy in the 
assumption concerning AV (i.e., AVi = AV 2 ) rather than the true difference 
in the time histories of the profiles. The time histories will be shown in 
later figures and the match will appear much closer. 

Because the required lift is within the capability of the configuration 
and the small difference in the acceleration profile, the approximate match 
appears attractive and is studied to the exclusion of the other approach 
during the remainder of this report. 

Computing the roll -angle profile of the flight vehicle requires the selec- 
tion of an appropriate set of initial conditions for altitude, velocity, and 
flight -path angle, in addition to specifying the flight vehicle L/D. In the 
next section, we will discuss the manner in which these quantities are 
selected so that the maximum acceleration and maximum acceleration rate 
encountered during the initial dive into the atmosphere by the flight vehicle 
will match the respective value of the reference vehicle. It will then be 
shown that if the flight vehicle has the aerodynamic lift capability to match 
the initial dive portion of the acceleration profile of the planetary mission 
vehicle entry, then a match of the subsequent constant altitude deceleration 
portion is also within its capability. 


Determination of Initial Conditions 

A vehicle entering the earth* s atmosphere will experience an acceleration 
buildup and a peak acceleration directly related to the initial velocity, 
flight -path angle, and vehicle lift-to-drag ratio. It can be shown that for 
a specified entry velocity and vehicle L/D, a unique correspondence exists 
between the peak acceleration rate and the peak acceleration encountered 
during the entries throughout a range of initial flight -path angles. This 
relationship between At^x and At^g^ as a function of L/D is shown in 
figures 4(a) and 4(b) for two values of Vq. These data were generated on 
the apalog computer simulation of the initial dive. The peak values of Aq 
and A-q were obtained during entries initiated at atmospheric encounter (the 
altitude at which 0.05 g is encountered) for a range of initial flight -path 
angles from -3° to -11°. Figure 4(a) presents these data for entries from 
earth orbit (Vq - l); figure_4( b ) shows the same type of data for entries at 
twice the orbital velocity (Vq = 2.0). 
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(a) Return from orbital mission. (b) Return from planetary mission. 


Figure 4.- Acceleration and rate of change of acceleration attained during earth entry. 

For entries from earth orbit at constant L/D several acceleration peaks 
may occur during the course of the entry. Only the peak At and occurring 
during the first acceleration buildup are plotted in figure 4(a). Our inter- 
est is in matching the peak. and A^_ for the two_ types of entries. In 
certain regions A t jrnQV and -E+ for entries at V-? = 1 and 2 can be related 

through the variable L/D and the initial flight -path angle. Two entries, 
one at = 1 (fig._ 4 (a)) and th£ other at - 2 (fig. 4 (b)), will achieve 

the identical peak and peak A^_ if L/D and 7 ^ bear the relationship 

shown in figures 4(a) and 4(b) . For example, to achieve an entry with a peak 
acceleration of 10 g and a peak acceleration rate of 0.5 g/sec requires that 
a vehicle entering at V n * = 2 have L/D - 0.2 and 7 - -5*8° (from fig. 4(b)). 

A vehicle entering at Vi - 1 would require that L/D =0.5 and 7 = - 7 .O 0 
(from fig. 4(a)). For this case, the orbital entry vehicle requires a greater 
aerodynamic lift capability than th§ planetary mission entry vehicle to 
achieve the identical peak At and At . For other cases, the L/D required 
may be less for the orbital entry vehicle than for the planetary mission 
entry vehicle. An. example of this is an entry that achieves a peak At of 
4-5 g and a peak At of 0.12 g/sec. A vehicle_entering at Vq = 2 requires 
that L/D = -0-5, while a vehicle entering at V± = 1 requires that 
L/D = +0.35* The relationship between the lift requirements for a planetary 
mission type entry and. the corresponding orbital mission entry that attains 
the same peak At and At can 1 ° e seen directly by cross -plotting the 
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(L/D)x and (l/d) 2 values of figures ij-(a) 
and Ij-(b) for equivalent A tma ^, A tmx 
pairs- The resulting variations are pre- 
sented in figure 5* Each point on fig- 
ure 5 represents two entries, one for the 
orbital mission return (V^ = l) and the 
other for the planetary mission return 
(Vi = 2), both of which attain the identi 
cal peak acceleration and peak accelera- 
tion rate- The aerodynamic lift required 
during each entry is designated (l/d) x 
and (l/D) 2 , respectively. Thus, from fig 
ure 5, the appropriate L/D can be estab 
lished for an earth orbital mission 
vehicle sojthat it is possible to match 
the peak A-^ and A-j- encountered during 
the initial dive portion of the planetary 
mission entry vehicle entry. The trian- 
gular shaded region on this plot repre- 
sents entries in which the planetary 
mission entry vehicle does not enter the 
atmosphere sufficiently to be captured, 
but skips back out. Points along the 
border of this shaded region represent 
entries at the overshoot boundary. An 
overshoot boundary exists for entries at 
+L/D (lift vector up) as well as the more 
commonly understood negative overshoot 
boundary ( -L/D, lift vector down). For 
entries at the positive lift overshoot boundary, the problems for vehicle 
control are compounded; for, in addition to skip out, there is an additional 
problem of exceeding the maximum allowable peak g. The vehicle enters with 
lift vector up (cp = 0 ) to keep the peak g within bounds. At pull out, a 
roll-over maneuver directs the lift vector downward (ep = l80°) to attain cap- 
ture. If the lift force is insufficient, the vehicle will skip back out of 
the atmosphere. For each entry of this type, there is a minimum allowable 
L/D to insure capture- These values of L/D and the corresponding peak 
acceleration define the +L/D overshoot boundary. For entry at the negative 
overshoot boundary, the peak accelerations will be low and of little concern. 
The vehicle enters with lift vector down (cp = l80°) and retains this attitude 
through the pull out. 



Figure 5.- L/D requirements for the initial 
dive phase. 


The undershoot boundary appears in figure 5 at the l6 g line. The under- 
shoot boundary, which specifies the maximum allowable acceleration, was set 
arbitrarily at l6 g for this study, but could have been set at any reasonable 
value. A 10 -g boundary is more commonly used. Two entry examples given 

earlier, one with A^ - 10 g, the other with Af = 4.5 g ; are repre- 

~max max 

sented in figure 5 by the square and triangular symbols, respectively. It is 
seen that the two examples actually represent entries at the positive and 
negative overshoot boundaries. An additional example at the l6-g undershoot 
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boundary is also indicated in figure 5 by the circular symbol. The figure 
shows that the acceleration profile of a planetary mission entry vehicle 
(L/D -0.4) entering at the l 6 -g undershoot boundary may be duplicated by an 
earth orbital mission entry vehicle with L/D = 0.5* The appropriate initial 
flight -path angles for these entries are 72 & -6.7 and 71 ^ -9-4. The values 
of 7 ^ are determined from figures 4(a) and 4(b) . 


Results of the Matching Process 

The three examples indicated in figure 5 will be used to illustrate the 
results of the matching process. The technique for matching the acceleration 
profiles was explained in the analog computer section of this report. The 
results of these entry simulations appear in figure 6, where the entries are 
shown from the time of atmospheric encounter until the flight vehicle reached 
an altitude of about 14.3 km (47,000 ft, 3h = 2) . This period corresponds to 
the initial dive phase and a portion of the constant altitude deceleration 
phase of the reference trajectory. The altitude, control, and acceleration 
profiles are shown as a function of time for the reference and flight vehicle. 
Zero time is chosen as the time at which the flight vehicle encounters an 
acceleration equal to 0.05 g* 

In each figure, the altitude profile shows the flight vehicle diving ini- 
tially to a lower, more dense altitude and remaining below the reference tra- 
jectory altitude throughout the entry. During the terminal phase, the flight 
vehicle dives progressively steeper into the atmosphere while the reference 
trajectory is a constant altitude. 

The acceleration buildup is slightly more gradual for the flight vehicle 
than for the reference, which is characteristic of a slower vehicle. The 
magnitude and time of the acceleration peaks coincide. The acceleration rate 
peaks occur at approximately the same time and are of the same magnitude. 

(The point of occurrence of the peak acceleration rate is not readily discern- 
ible on these acceleration plots.) Good agreement between the acceleration 
profiles is achieved up to the point where the flight vehicle can no longer 
generate sufficient aerodynamic lift to pull down into the atmosphere at the 
progressively increasing rate required. Consequently, the flight vehicle 
acceleration profile diverges from the reference acceleration profile after 
this point, at about 94 seconds for the l6-g entry, l44 seconds for the 10 -g 
entry, and 326 seconds for the 4.5-g entry. 

The roll control profiles are shown in the lower portion of each figure. 
For the flight vehicle in each entry, L/D = 0-50. The L/D for the reference 
vehicle is 0.4, 0.2, and 0.5 for the l6-g, 10-g, and 4.5“g entry, respectively. 
In addition to the actual control profiles for the reference and flight vehi- 
cles, the nominal roll control profile for the flight vehicle is shown. These 
nominal roll control profiles (as explained before) are the quantitative 
solutions of equation ( 13 ) when the_values of the planetary mission entry or 
reference vehicle state variables (A^, 7 2 , and (l/d) 2 ) are used as the refer- 
ence values. The control system of the flight vehicle attempts to control 
about this nominal. Excursions from the nominal are due, in part, to the 
closed-loop nature of the control and, in part, to the constraints on the 
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control action. One constraint limits the vehicle roll rate to less than 
±20° per second. This caused a slight deviation of the control profile at 
roll over and during the terminal phase of the simulation where the nominal 
specifies a higher roll rate. An additional constraint does not allow the dip 
or discontinuity that appears in the nominal profile just prior to roll over 
in figures 6(a) and 6(b) . This results in a smooth and continuous roll-over 
maneuver similar to that of the reference vehicle control profile instead of 
the discontinuous one specified by the nominal. The effect of this constraint 
on the acceleration profile is not noticeable. 

Figure 6(c) shows the example time history for entry at the negative over- 
shoot boundary. This case merits further discussion because it differs from 
the other two examples for entry at the positive overshoot and undershoot 
boundaries. During the initial dive, the flight vehicle enters with a roll 
angle of about 42°, indicating that the full lift capability of the vehicle is 
not required to achieve the match. This was not the situation for the other 
two examples which both required that the lift vector be full up (cp => 0°) 
during the initial dive. It can be shown to be generally true that the (l/D) x 
required to simulate the entry at the negative overshoot boundary is less 
than the (l/d) 2 required for the actual planetary mission return entry. This 
can be verified from figure 5 for the example given here for At max - 4.5 By 
as well as for many other examples that fall along the negative overshoot 
boundary. The opposite result (i.e., (l/D) x > (l/D) 2 ) is generally true at 
the positive overshoot and the undershoot boundaries. 

In figure 6(c) the shape of the 
flight vehicle roll control profile 
for entry at the simulated negative 
overshoot boundary does not resemble 
the control profile for the corre- 
sponding reference vehicle entry. 
There is, however, a close similar- 
ity in the shape of the reference 
and flight control profile for the 
other two cases (figs. 6(a) and 
6(b)). 


Simulation Time 


The maximum time duration over 
which a simulated planetary mission 
entry may be sustained is presented 
in figure 7* Lines of constant time 
are given as a function of the peak 
acceleration encountered during the 
entry and the vehicle L/D. The 
values represent the total time 
from atmospheric encounter 
(A-f- =• 0.05 g) until the flight vehi- 
cle can no longer sustain the 
desired level of acceleration. The 
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Figure 7.- Maximum duration of planetary entry 
simulation; time - 0 at - 0.05 g* 
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regions corresponding to the planetary mission return overshoot and undershoot 
trajectories are shown in this figure as was done in figure 5 - The points 
corresponding to the three entry examples given previously (see figs. 5 and 6 ) 
are also included for reference. 

It is apparent from the shallow slope of the curves that the time duration 
of the simulation is not a strong function of the vehicle L/D. For instance., 
a vehicle with an L/D of 0-50 could sustain a simulated planetary entry for 
l44 seconds for the = 10 g case. This is indicated by the square sym- 

bol on figure J . An increase in the vehicle L/D to O. 7 O (40-percent 
increase) would result in an additional 20 seconds (l4-percent increase) of 
simulation time for this 10 -g case and would not affect the ability to attain 
a match during the high acceleration portion of the entry. 


Aerodynamic Heating 

For the entries from earth orbit considered in this study, only the aero- 
dynamic heating that results when a laminar boundary layer is assumed will be 
evaluated. The Reynolds number can be shown to be less than 3xl0 5 at peak 
heating for each of the entries of interest, making it reasonable to expect 
a considerable extent of laminar flow. The laminar convective heating rate 
per unit area at the stagnation point is given by 

q = Kv/p V 3 

The constant K is selected so that the heating rate at peak heating matched 
that given by a presumably more accurate method of computation for the entry 
of the same configuration. In this study , the entry trajectory of the Apollo 
204A mission was simulated on the analog computer using conditions as given in 
reference 5- The constant K was selected to match the peak heating rates 
for this entry. Results of the computations of the maximum heat rates and the 

total heat load defined as 

H = I q dt 
°o 

are presented in figure 8 and com- 
pared with data obtained from ref- 
erence 5- The orbital mission entry 
trajectories of interest in this 
study fall in a region near the 
origin in this plot. The three entry 
entry examples used previously are 
indicated within the region by the 
circle, square, and triangle for 
l 6 g, 10 g, and 4.5 g? respectively. 
The peak heating rates and total 
heat loads experienced are much less 
than the capability of the Apollo 
capsule (lunar return configuration) . 
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Figure 8.- Aerodynamic heating at the stagnation 
point . 
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The total heat loads are also less than the typical Apollo orbital flight 20 4A, 
although the heating rates are higher. Data points for Apollo flights SA201 
and SA202 are also presented for comparison. These two flights are for 
ballistic lobs at higher entry velocities (28,500 ft/sec). 


Range Control Considerations 

During the entry from orbit, a spacecraft would perform range control 
maneuvers in order to arrive at a desired destination. Ideally, range error 
would be corrected most efficiently during the early portion of the entry 
where the maneuver capability of the vehicle is greatest. Performing range 
control during the early portion could conflict with the procedures required 
to fly the acceleration profile that was desired for the planetary entry simu- 
lation. To avoid this complication, control to achieve the desired accelera- 
tion profile and control to achieve a desired landing point are separated into 
two distinct phases in this study. A resulting decrease in range control 
capability will be one of the costs in performing this type of entry. Con- 
sider then that at some point during the entry, the planetary entry simulation 
will be terminated by the earth orbital mission return vehicle and terminal 
range control initiated. We will direct our attention to the effect simula- 
tion duration has on the range capability of the entry vehicle as specifically 
indicated by the size of the landing footprint . Using the previous example of 
the entry at the positive overshoot boundary that results in a maximum accel- 
eration of 10 g, figure 9 shows the situation in which the control maneuver 
required to match the acceleration profile is terminated at a time prior to 

the absolute end point for the simula- 
tion. Three trajectories are shown in 
figure 9, two of which illustrate 
range control; the third trajectory is 
identical to the trajectory presented 
in figure 6(b) . The trajectory desig- 
nated Q) indicates a situation in 
which the simulation is terminated 
after 80 seconds by rolling the vehi- 
cle so that the lift vector is up 
(cp - 0°) • This maneuver reduces the 
total acceleration on the vehicle as 
it rises in the atmosphere, and, ulti- 
mately, the vehicle is carried farther 
downrange. The second example (desig- 
nated (2) ) indicates that the vehicle 
has rolled so that the lift vector is 
down. In this case, the vehicle dives 
deeper into the atmosphere, causing an 
increase in the acceleration and ulti- 
mately shortening the range to touch- 
down. During the analog computer 
analysis this procedure is carried out 
not only for the full -up and full -down 
direction of the lift vector, but for 
all directions between ±l80° of roll 
angle. The trajectory that results for 



Figure 9*- Entry trajectory characteristics with 
range control; A raax = 10 g, L/D = 0-5- 
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shorn in figure 10. 


each of the constant roll angle conditions is computed down to the terminal 
points (at 14-3 km , Ph = 2), which are plotted to form a landing footprint as 

The footprint designated 80 seconds corresponds to the 

situation just described in which the 
simulation is terminated 80 seconds 
after entry. Progressively larger foot 
prints would result if the simulation 
were terminated at earlier times , with 
the largest area representing the situ- 
ation where there is no simulation at 
all (t — 0 ) , and the ground area attain 
able is that resulting from using the 
maximum lift capability of the vehicle 
from the initial point at atmospheric 
encounter (At =• 0.05 g in these exam- 
ples) down to the terminal point of 
14. 3 km altitude for each case. 

For each trajectory in figure 10 
the acceleration peak was at least 10 g 
and, for some of the trajectories in 
which the roll angle is greater than 
90 °, a second acceleration peak greater 
than 10 g (see @ on fig. 9) but 
always less than 20 g occurs . 



Figure 10.- Attainable ground area after 

termination of p lane tar y_entry simulation 
of different durations; A^y = 10 g, 

l/d = 0.5. 



Duration of simulation, sec 

Figure 11. - Vehicle maneuver capability as a 
function of elapsed time ffom atmospheric 
encounter (At = 0.05 g) - 


It is apparent in figure 10 that 
there is a rapid decrease in footprint 
size or vehicle maneuver capability for 
simulation durations in excess of 
50 seconds. This is readily apparent 
if the area of each footprint in fig- 
ure 10 is determined and plotted as a 
function of simulation duration. In 
figure 11, these data are shown in 
dimensionless form for the 10 -g case 
discussed here as well as for the cases 
that correspond to a 4.5“g entry and a 
l6-g entry (discussed previously). The 
high cost, in terms of loss in landing 
area available, for increased simula- 
tion time is indicated by the three 
curves in figure 11. 

Midcorridor Entries 

To illustrate the capability of an 
earth orbital mission vehicle to match 
the acceleration profile corresponding 
to the possible midcorridor entries of 
a planetary mission return vehicle, a 
summary plot similar to that in 
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figure 5 of the control requirements 
for the initial dive portion is pre- 
sented in figure 12. Figure 5 pre- 
sented the capability for simulating 
entries at the overshoot and undershoot 
boundaries for different values of L/D, 
while the curves of figure 12 are for 
entries in the midcorridor range for 
only one value of L/D. Each point on 
the plot represents two entry trajecto- 
ries, one for Vi = 2, the other for 
= 1. Both trajectories attain the 
identical peak acceleration and peak 
rate of change of acceleration during 
the initial dive phase. Whereas 
entries at the undershoot and overshoot 
boundaries require the lift vector of 
the planetary mission entry vehicle to 
be either full up (cp = 0°) or full down 
(cp = l 80 °), the midcorridor type 
entries allow a component of the lift 
vector in the lateral plane. The vehi- 
cle roll attitude during the initial 
dive would still be constant but at a 
selected value between zero and ±180 . 

To obtain the data of figure 12, the flight and reference vehicle are 
both ' considered to have an L/D of 0.5* A slightly different plot would 
appear for each combination of and (l/D) 2 . Figure 12 is for just one 

of these combinations . 

The midcorridor entries on figure 12 are represented in terms of the peak 
acceleration attained and the roll attitude maintained by the reference vehi- 
cle, cp 2 , and required of the flight vehicle, cpi . The required flight -path 
angle for each vehicle has been omitted for the sake of clarity even though it 
is important in establishing the initial conditions. 

As an example of the simulation of a midcorridor entry, consider a plane- 
tary mission vehicle entering the earth atmosphere at a roll attitude of 90 ° 
and at the appropriate flight -path angle for attaining a peak acceleration of 
8 g. The peak acceleration and peak acceleration rate could be matched by an 
earth orbital mission vehicle entering at the appropriate flight -path angle 
with cpi =■ 35°* 

The example given previously corresponding to entry at the negative over- 
shoot boundary is also indicated on figure 12 (triangle symbol) . This is the 
entry example illustrated in figure 6(c), which attains a peak acceleration of 

4.5 g. 

Numerous other examples for simulation of midcorridor -type entries are 
within the capability of the orbital entry vehicle for the example given in 


^•mox * *3 



Figure 12.- Control requirements during initial 
dive for simulating raidcorridor-type entries; 
(L/D)i = 0.5, (L/D) 2 = 0.5. 
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figure 12 [(l/d) x =0-5, (l/D) 2 - 0.5]- All entries that fall between the two 
shaded regions on the plot represent feasible midcorridor entry simulations. 


Retro Velocity Requirements 


Thus far, no consideration has been given to how the entry vehicle 
arrived at the edge of the earth atmosphere flying at the correct flight -path 
angle to give the acceleration buildup and peak g anticipated. Under the 
assumptions that have been made for the initial dive into the atmosphere 
(Vi = 1 and cp = constant), the peak g and maximum acceleration rate are a 
function of. the L/D of the vehicle and the initial flight -path angle only. 
The proper L/D is determined as explained previously. The flight -path angle 
at entry will be a function of the initial orbit parameters and the retro 
velocity increment at the time of the deorbit maneuver. 




.04 .08 

Retro velocity, SV 


Figure 13.- Retro velocity requirements. 


The curves presented in figure 13 
show the relationship between the 
flight -path angle that would be real- 
ized at entry and the retro velocity 
increment added during the deorbit 
maneuver for three different initial 
orbital altitudes . On this plot it is 
possible to show a region of interest 
corresponding to planetary mission 
entry simulation and to compare this 
region with the available retro capa- 
bility of the Apollo orbital vehicle. 
The vehicle is considered to be in a 
circular orbit prior to retro. The 
Apollo orbital retro velocity increment is added in a 
^configuration direction to give, the maximum negative 
flight -path angle when the vehicle 
arrives at the atmosphere. 


A second ordinate axis is included 
on figure 13 to show the approximate 
peak acceleration that the vehicle 
would encounter upon entering with the 
lift vector up (cp = 0°).. The values 
At^ax would be the lowest peak acceler- 
ation encountered for a designated y±. 


Higher peak accelerations could be encountered for entry at a roll angle 
greater than zero. The points corresponding to the three example entries are 
again included on this plot for trajectories from an initial orbital altitude 
of 322 km (200 statute miles). 


It may be noted that for the 4-5 g example (triangle symbol) the vehicle 
does not enter with cp = 0° (see fig. 6(c)), and thus this data point does not 
relate to the acceleration ordinate on figure 13 - For the other two examples, 
l6 g and 10 g, the roll angle is zero upon entry and the values given on the 
acceleration ordinate are correct. 
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The region of interest for planetary mission entry simulation is shown on 
the plot. Entries that result in peaks between 10 and 1 6 g, simulating 
entries at the undershoot boundary, show that large retro velocity increments 
are required in the range from 6V = 0.12 to 5V = 0.17 (950 to 13^-0 m/sec) for 
return from a 322 -km orbit. The entries for midcorridor and overshoot bound- 
ary simulation (cp > 0°) require smaller y± or 5V and would extend the lower 
limit of the region of interest to include that shown in the figure (&V - 0.01 
to 6V = 0.17). It is noted that part of the region of interest lies beyond 
the capability of the present Apollo orbital configuration. 

CONCLUSIONS 


(1) The initial acceleration buildup (up to peak g) experienced during 
atmospheric entry of a planetary mission return vehicle may be duplicated dur- 
ing the entry of a manned earth orbital mission through proper choice of the 
lifting capability of the entering spacecraft and establishment of the appro- 
priate initial flight -path angle. Accurate duplication of both the maximum 
acceleration and the maximum acceleration rate can be achieved for entries 
that simulate planetary return at the undershoot boundary, positive and nega- 
tive lift overshoot boundaries, as well as midcorridor approach with a 
capsule -type spacecraft configuration. 

(2) The acceleration profile corresponding to the constant altitude por- 
tion (after peak g) of a planetary mission return entry may be duplicated, 
in part, through proper control of the direction of the lift vector of the 
entering orbital spacecraft . 

( 3 ) The lift requirements for the orbital return vehicle are generally 
determined by the requirement to achieve a match of the acceleration profile 
during the initial dive into the atmosphere prior to peak g. The lift 
required to simulate the constant altitude portion of a planetary entry is 
generally less than this initial segment. Only small and inefficient gains 
result from increasing the vehicle lift capability above that required during 
the initial dive in order to increase the duration of the simulation of the 
constant altitude segment . 

(4) Dowrrrange and crossrange capabilities are reduced by control maneu- 
vers required to achieve the acceleration profile match. The subsequent 
reduction in the size of the vehicle landing footprint is a strong function of 
the velocity (or time) at which the simulation is terminated and range control 
initiated. 

(5) Large entry flight -path angles, implying large retro velocity incre- 
ments, are required by the orbital -type entry vehicle to simulate a high g 
planetary mission entry at the undershoot boundary. Lower retro velocity 
increments are required for simulating midcorridor and overshoot boundary type 
entries . 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, June l4, 1968 
125-19-01-42-00-21 
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APPENDIX A 


CONDITIONS FOR USING EQUATION ( 8 ) 


To obtain a quantitative solution of equation (13) in the analytic 
section, we have considered that the total velocity change experienced by the 
reference and flight vehicles during the entries was identical at correspond- 
ing times along the trajectories and, further, that this total velocity change 
was equal to the integral of the total acceleration vector: 




dt 


or AV 2 - AVi 


(Al) 



In sketch (a) the change in total velocity 
for each vehicle is 


AV = — + g sin 7 At 

v m ° ' 


(A2) 


or, in dimensionless integral form, 
' T ' "At 


AV = 


s/TT (l/d ) 2 


+ sin 7 ) dT 1 (A3) 


where the drag acceleration has been replaced by the equivalent value in terms 
of the total acceleration = (D/mg) \l 1 + (L/D) 2 . 

From equation (A3) in order for the velocity change to be equal for the 
two vehicles, a unique correspondence must exist between the flight parameters 
A-fc, L/D, and 7 for the two vehicles. Since the flight vehicle _is controlled 
so that the acceleration A^ is equal to the reference value A-fc^ at each 
point along the trajectory, it is apparent that an error would develop if an 
inequality exists in the values of (L/D)i, (l/D) 2 , 7 ±> and 72 * This error 
would be 


+ (sin 7 2 - sin 7 i)dT 

(AU) 

Although it is feasible to have (l/D)i =• (l/d) 2 , it should be clear that 
the flight -path angles could not be equivalent at each point along the trajec- 
tory (see altitude time histories, fig. 6 ). Thus an error will always develop 
as a result of this difference in flight_-path angles. Depending on the partic- 
ular trajectory flown, this error in A V may vary from 2 to 9 percent of the 
total velocity change at the end of simulation. This error will have only a 
small effect on the acceleration matching results. 
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APPENDIX B 


TIME EQUATION 


An approximate equation may be derived that relates the accumulated time 
along the entry trajectory with the velocity change for the constant altitude 
portion of the entry. For the constant altitude phase of the reference vehi- 
cle entry (7 0° and p = const) , equation (A2) reduces to 


AV - 



(Bl) 


from which the change in time 
ity as 

t - t j 


can be related directly to the change in veloc- 



dV 


Vj At/ n/TT (L/D ) 2 


(B2) 


Time 



where the limits of integration extend from the 
velocity at initiation of the constant altitude 
(pull out) to the point of interest along the 
trajectory (designated * in sketch (b)). 

For this constant altitude phase (constant 
spheric density) , the total acceleration varies 
the square of the velocity. 


At = KV 5 


where 



phase 


atmo- 

as 


(B3) 


Thus, we can substitute for the acceleration and obtain an expression that can 
be integrated: 


t - ti - £ Ai-LiiMid f v av 

K JtT V 2 


Vi v 

After integrating and substituting the limits, we obtain 


(Bit) 


t = - 


s/i + (l/d) 2 2 /l 1 


K 


y " v T J + 1]: 


(B5) 


or , in terms of the independent variable AV for an initial velocity of 
V 2 = 2, 




sFl + (l/d) 2 2 


k 


- AV 


AVj. 


W~ + t I 


(b6) 
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APPENDIX C 


EQUATIONS OF MOTION 


Equations for the angular and translational motions of the spacecraft 
(three translational degrees of freedom and one rotational degree about the 
roll axis) are given herein. The equations are given in terms of a modified 
flight -path axis system (H frame) that is rigidly attached to the vehicle and 
moves with it. The three translational equations (from ref. 3) are 

u h w h x h 


— (motion is in equatorial plane) 


Wh + -f - So 


where Uh and Wh are the vehicle velocity components in the H frame and 
Vh =■ 0 by definition of the H frame (see ref. 3) • 

To simplify the equations with respect to computer scaling, the following 
dimensionless variables are introduced: 


Ap = — - 1 
^ r 0 

Uh - u hr 


(altitude) 


(velocity) 


T= ^t 


(time) 


In terms of these new variables, the translational equations of motion become 


1 + Ap 


' T X h 

(l + Ap) — ^ dT 1 - Ap 

, mg D 


1 + Auh mg 0 


Ap - 
A + Ap) 2 


2Au h + (Au h )‘ 
1 + Ap 


h _ i ■■ i ■■■■■■ ■■■■!■ iii ini ■ ■■iimi«iii ih in w miw i an ■■ ii mil inau 


n«i ■ 1 1 1 1 i I II I 


The external forces on the vehicle (not including gravity) are 


x h D / L . 

= - — sm y + cos y 

mgo mg 0 \ D 


Y h L f x h 


tng 0 D \rago. 


sin cp 


z h D (h ^ . 

= — cos y + sin y 

“So “So VP 


where the flight -path angle 7 is defined as 


7 =■ tan" 


-1 r «h 


The total acceleration and rate of change of acceleration is 


A t = — 1 + ( | Y 

15 “So s vVmax 


^ ( 2 D/“Sq 

t ^ -Vug + w| 


Auxiliary relationships are: 


Total velocity = (U§ + W § ) 1/2 


Downrange 


= ^ (l + Aufr)cos ^ dT* 


Crossrange = / (l + Au^)sin ^ dr* 
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